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Within the geometric optics approximation, we derive the 
light cone condition for a class of homogeneous non-trivial 
QED vacua using the effective action approach. Our result 
generalizes the "unified formula" suggested by Latorre, Pas- 
cual and Tarrach which turns out to be a low-energy limit. We 
calculate the velocity shifts induced by electromagnetic fields 
and finite temperature and especially analyze the high-energy 
limit consistent with the loop approximation. 
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I. INTRODUCTION 

The vacuum of classical electrodynamics is denned by 
the absence of charged matter: j M = 0. Since the the- 
ory is linear, this vacuum is unique and trivial. As an 
immediate consequence, the light cone condition which 
characterizes the propagation of light also is uniquely de- 
termined and trivial: k 2 = 0. 

Considering the (more) fundamental theory of quan- 
tum electrodynamics, we can only demand the absence 
of external currents j£ = 0. But non-charge-like modi- 
fications of the vacuum such as external fields, gravita- 
tion, temperature or non-trivial topology can influence its 
ubiquitous quantum fluctuations by acting on the prop- 
erties of the fluctuating fields (charges, masses,. . . ). In 
general, the expectation value of the electromagnetic cur- 
rent will hence not vanish if it is taken with respect to a 
vacuum that is exposed to a modification z. According 
to Schwinger Q , we can calculate the induced current to 
one-loop order with the aid of the formula 



J' 



-ie tr 7 M G(a;, x; z) 



(1) 



where G(x,x;z) denotes the Green's function (electron 
propagator) in presence of the modification z evaluated 
at the same space-time point x (closed loop). Since prop- 
agating light may couple to this current, the light cone 
condition alters k 2 7^ and the vacuum is called non- 
trivial. Light velocity shifts of this kind were found by, 
e.g., Adler ||, Brezin and Itzykson || for magnetic fields, 
by Drummond and Hathrell [|J for gravitation, and by 



Scharnhorst [|| and Barton Q for a Casimir configura- 
tion. Further important examples are studied in refs. 

In this letter, we investigate a general approach to ob- 
tain the light cone condition for a class of non-trivial 
QED vacua. A first unifying result has been given by 
Latorre, Pascual and Tarrach ||, who identified the so- 
called "unified formula" for the velocity shift 
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(2) 



where m denotes the electron mass, u the (renormal- 
ized) background energy density and a ~ 1/137. Indeed, 
eq.(0) perfectly describes the polarization and direction 
averaged findings of refs. [g [K]] in the low-energy do- 
main. (For a gravitational background, one a has to be 
replaced by the combination (G^m 2 ) involving Newton's 
constant). In the case of gravitation, Shore (T^| proved 
a polarization sum rule that represents a generalization 
of eq. (||) . Additionally, he studied the case of weak elec- 
tromagnetic fields in a consistently covariant manner (cf. 
the electromagnetic birefrigent part of ref. j7|). 

In the present work, we address the question of whether 
a generalization to arbitrary vacua is possible. We es- 
pecially aim at the inclusion of the high-energy domain 
where eq.(||) obviously fails, since the energy density is 
naturally unbounded from above. 



II. LIGHT CONE CONDITION 

Instead of calculating the current expectation value 
eq.(0), we make use of the effective action approach. 
Once the high-energy degrees of freedom of the full quan- 
tum theory are integrated out, the effective action by 
definition describes the non-trivial vacuum as a classical 
medium. The following assumptions are essential for the 
formalism and classify the type of vacuum under consid- 
eration: 

1) The vacuum modification is homogeneous in space and 
time (or at least slowly varying compared to the wave- 
length A of the propagating light f 1 ")- 

2) The wavelength of the propagating light is large com- 
pared to the Compton wavelength (soft photon approxi- 
mation) : Lu/m <C 1 . 

3) Vacuum modifications caused by the propagating light 
/ M " itself are negligible. 
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4) Vacuum modifications EM fields) behave passively 
towards EM fields. 

Referring to assumption 1) and 2), we can neglect any 
derivative term of the field strength in the effective La- 
grangian. In particular, assumption 1) is equivalent to 
the geometric optics approximation and assumption 2) 
excludes dispersive effects from the formalism. 

Assumption 3) justifies a linearization of the field equa- 
tions with respect to /'*". 

The meaning and necessity of assumption 4) will be 
explained when we resort to it. Note that we demand 
neither a low-energy modification of the vacuum nor that 
the deviation from the Maxwell action be small. Besides, 
it is understood that we only take into account the real 
part of the effective action, assuming that the vacuum is 
sufficiently stable. 

First, we consider a vacuum that is purely modified 
by an electromagnetic background field. Restricted by 
assumptions 1) and 2), the effective Lagrangian can only 
depend on the two Lorentz and gauge invariants of the 
Maxwell field 

C = C(x,y), (3) 

where we introduced the linearly independent invariants 

x:=\F^ v = \{B 2 -E 2 ) (4a) 
y := \F MV *F» V = E-B. (4b) 

The field strength and its dual are defined as usual: 



= l e ^ af3 F al3 . 
We obtain the field equations from C by variation 
dC 



(5a) 
(5b) 



= 8. 



d ll (d x £F» t/ + d y £*F^) 

F a0 , (6) 



where d x ,d y denote the partial derivatives with respect 
to the field strength invariants (0) and M£g is given by 



<; := F^F a0 [dlL) 



*F^*F a0 (d 2 C) 



-d xy C (F» v *F a g + *F flv F a ^ 



(7) 



To arrive at the desired light cone condition in the spirit 
of Shore's covariant formalism we make use of the 
following five pieces of information: 

(i) We split up the field strength into its constant (or 
slowly varying) background part and the propagating 
part F^ -> F^ v + f^, and linearize with respect to f» v 
(assumption 3)). In Fourier space, the field derivatives 
are represented by d^F KX — > fc M / KA , where denotes 
the wave vector of the propagating light. 

(ii) In Fourier space, f^ v can be written in terms of the 
polarization vector e^: f 1 " cx (k^t v - fc'V). Without 
loss of generality, we choose the Lorentz gauge fc M e M = 0. 



(iii) The average over polarization states can easily be 
taken with the aid of the well-known rule X) po z e ^ eU ~ * 
gP", where we use the metric g = (— , +, +, +). Note that 
the additional terms on the RHS vanish by the antisym- 
metry properties of M^ g . 

Intermediately, the field equation @ yields 



= 2{d x C) k 2 + MZ k^k a 



(8) 



(iv) Using the fundamental algebraic relations of the field 
strength tensor and its dual (T) 



F tJ,a F u a ~ *F^ a *F v a = 2 x tT 
F^ a *F v a = *F^ ia F u a = y gT , 



(9a) 
(9b) 



the Lorentz structure of can be decomposed into one 
term proportional to 5£ and another proportional to the 
Maxwell energy-momentum tensor 



rpll 
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(10) 



(v) Only the vacuum expectation value of the energy- 
momentum tensor is physically meaningful and can be 
defined by a variation with respect to the metric 



-T^{d x C) + g^ (C - xd x C - yd y C) . (11) 



Solving eq. ( |ll| ) for T^ u and inserting into eq. (|10|) , we can 
present in its final shape 



M ^ = 2 



' l (d 2 +d 2 )C 
2 d x C 

+yd xy C+ ~ 



(T^) xy + S^ 2 x(d 2 -d 2 y )C 



2 ' d 2 v )c 

(L-xd x £-ydy£) 



d x C 



(12) 



Substituting into eq.(g), we arrive at the desired 

light cone condition for EM field-modified vacua fulfilling 
the above-mentioned assumptions. 



where 



Q = 



Q {T^ ) X y h^hjy , 

\{dl + dl)C 



(13) 



(9*C) 2 +(9«£)(§ (dl-dl)+yd^C+%dl+dl)C(l-xd*-yd y )C 



(14) 

Note that the Lorentz structure of the light cone con- 
dition eq.(|l3|) is identical to the one found by Shore in 
the limit of weak electromagnetic fields JTl[| . Hence, it is 
the Q-factor that will essentially govern the high-energy 
domain. 

Regarding the universality of the "unified formula" 
eq.(2], we are aiming at an extension of the validity of 
eq. ( 13) to arbitrary non-trivial vacua. Therefore we take 
the expectation value of the light-cone deforming quanti- 
ties on the RHS of eq. (113) with respect to the additional 
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vacuum modifications parametrized by the (collective) 
label z 

k 2 = z (o\Q(Tn xy \o) z k^ 

= ^2 Z (0\Q \%% z {i\{T^) xy |0), Kk v , (15) 



condition is a generalization of the "unified formula" of 
Latorre, Pascual and Tarrach ||. However, the Q-factor 
generally depends on all the variables and parameters of 
C and we will have to verify that its low-energy value 
approaches the constant pre- factor of eq. (0) . 



where we inserted a complete set of intermediate states 
in the last line. Falling back on assumption 4), we make 
use of the fact that the vacuum should behave passively 
with respect to EM fields. In other words, once in the 
vacuum state of z, switching on EM fields should not alter 
this ground state: z (0| (x,y) \i) z = {(x,y)) z S i- (From a 
rigorous point of view, the application of our formalism 
to the case of gravitation drops out at this stage. In fact, 
we will only consider non-gravitational modifications in 
the following.) Since Q solely depends on x and y (via 
C(x,y)), this prohibition of a backreaction leads to 



Q\i)z = (Q) z 5 0t 



(16) 



Evaluating the expectation value of Q that is functionally 
dependent on C(x, y), leads back to the definition of £ via 
the functional integral over the fluctuating fields. E.g., if 
the modification z imposes boundary conditions on the 
fields, the functional integral has to be taken over the 
fields which fulfil these boundary conditions. This defines 
the new effective Lagrangian characterizing the complete 
non-trivial vacuum 



(Q) z = (Q(L(x,y))) z = Q(C(x,y;z)) 



(17) 



We finally end up with the light cone condition for a 
class of non-trivial vacua reconcilable with the above- 
mentioned assumptions 



k — Q{x, y, z) (T^ )xyz k^ku 



(18) 



As an exact statement in the sense of effective action 
theories, the validity of eq.(|l8|) is not restricted to per- 
turbation theory. 

Further useful representations of eq.(|l8|) are obtained 
by choosing a certain reference frame and introducing 



(19) 



Here we defined the phase velocity by v := fc°/|fe|. For 
eq.(|I^), we obtain 



v 2 = 1 - Q (T^kph 



(20) 



Averaging over propagation directions, the light cone 
condition yields 



v 2 = l-^Q{T™) = l-^Qu, 



(21) 



where we assumed that Q(T 00 ), (T a a ) < 1 and u denotes 
the (renormalizcd) energy density of the modified vac- 
uum. These representations indicate that the light cone 



III. APPLICATIONS TO THE LIGHT CONE 
CONDITION 

Up to now, our formalism has only worked in a sim- 
ple manner at the expense of dealing with effective La- 
grangians. Of course, the calculation of the latter gen- 
erally lacks simplicity. Indeed, we will only deal with 
effective Lagrangians of first-order perturbation theory, 
but the implementation of higher-loop corrections is di- 
rectly controlled by the effective action approach. Due to 
the perturbative properties of QED, first-order effective 
Lagrangians are appropriately characterized by 



C = C 



M 



c 



M 



«1, 



(22) 



where Cm, = —x denotes the Maxwell Lagrangian and 
C c contains the correction terms. For this class of La- 
grangians, the denominator of the Q-factor in eq.(|l4|) 
simplifies to 

denom.(Q) = 1 + 0{C C ) , (23) 
and the approximation Q — \ {d 2 + d 2 )C is justified. 



A. Weak EM Fields 

The weak field limit of the one-loop effective La- 
grangians of QED, i.e., the Heiscnberg-Euler Lagrangian, 
is given by 



C = —x + ci x 2 



where 



8a 2 



c 2 y 



Ua 2 



C-2 



(24) 



(25) 
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Employing eq. ( pl| ) , we immediately find the propagation 
and direction averaged velocity 



y = ci -+ 
v = l- 



C2 , 

44a 2 
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B 2 ) 



(26) 
(27) 



Here, we have rediscovered the "universal constant" of 
the "unified formula", eq.(||). Since all of the low-energy 
calculations are indeed based on the Heisenberg-Euler 
Lagrangian, the seeming universality of this factor is not 
astonishing. As will be pointed out in the following, 
this combination of constants only represents a natu- 
ral low-energy limit of the respective modifications. Of 
course, higher-loop corrections to the Heisenberg-Euler 
Lagrangian will also modify this factor. 



3 



B. Strong Magnetic Fields 

The vacuum modified by EM fields of arbitrary 
strength consistent with the one-loop approximation is 
described by Schwinger's famous Lagrangian JlJ 



100 

1 fds 



8n 2 j s- 
o 



~ ms (es) 2 \y\ coth(es(yJx 2 +y 2 + i)3 
cot\es(y x 2 +y 2 — x)^j — -(es) 2 x - 



(28) 

The convergence is implicitly insured by the prescription 
m 2 — ► m 2 — ie. For purely magnetic fields, i.e., on the 
positive x-axis in field space, we find the Q-factor by 
differentiation 



Q(h) 



IOC 

1 a f dz_ -2hz 

2a 2 7T J 

o 



z 



z coth z — 1 f 



sinh z 



-z coth z 



(29) 



where we introduced the convenient dimensionless pa- 
rameter /i := —: 4rS-. With some effort, the evalua- 
tion of the integral is analytically achievable. Details will 
be given in a forthcoming paper pjfl. Our findings are 



Q(h) = Ar- 



(2h 2 ~ + - ft - 3/i 2 - 4/ilnr(ft) 

(30) 



-2/iln27T 



where 'J denotes the logarithmic derivative of the F- 
function and £' is the first derivative of the Hurwitz Zeta 
function with respect to the first argument. 

The zero- field limit coincides with eq.(p6|). For strong 
fields, the last term of eq.(30) cx ^ oc |B|dominates the 
expression in the square brackets. Hence, the Q-factor 
decreases with 



Q(fl)~i^J-I,for 
6 7r B cr B 



B 



(31) 



Since the loop-corrections to (T^ u ) are of higher order 
and can be neglected, the velocity of propagating light 
will be modified according to 



f 



a sin 
" 2 



2h 2 



W(l+h)-4hlnT(h) - 3/i 2 



f 



2/iln27r+- + 4C'(-l,/i) 



1 

6h 
(32) 



One can show JL2| that eq.(|32|) coincides with the find- 
ings of Tsai and Erber jl(| . Although the velocity shift 
increases proportional to the magnetic field for large B 



(last term oc Jfc in eq. (J32J)) , its total amount remains 
comparably small 

o 

Sv ~ 9.58.. • 10~ 5 at B = B CI = — (33) 

e 

for strong B-fields consistent with the one-loop approxi- 
mation, i.e., -jp- < — ~ 430. In order to let Q(T' Al!/ )fc^fc i/ 
be bounded, we expect that higher loop corrections pro- 
mote the decrease of Q(B). 



C. Finite Temperature 

An effective Lagrangian describing a finite tempera- 
ture vacuum state can always be decomposed into its 
temperature-independent (T = 0) part and its thermal 
correction 



C(T) = C(T = 0) + A£(T) 



(34) 



whereby £(T = 0) denotes the usual zero-temperature 
Lagrangian, e.g., eq. (p8[) . In a similar manner, we can 
decompose the Q-factor 

Q(T) = Q(T = 0) + AQ(T) 

= \ [{dl + d 2 )£(T = 0) + {dl + d 2 ) A£(T)] . (35) 

Including a magnetic background field, A£(_B,T) was 
calculated by Dittrich 031 



A£(B,T) 



\Jii f ds 
4^ J ~J' 



s esB cot esB 



(36) 



xT 



e 2 (0,47risT 2 ) 



1 



2T\FWs 



where O2 denotes the second Jacobi O-function |H. 

Here, we encounter a subtlety of the formalism. The 
Q-factor is evaluated by deriving £ with respect to x and 
y. However, the finite-temperature formalism demands 
that the electric field vanishes in order to fulfil the prin- 
ciple of thermal equilibrium. Hence, we first have to pre- 
serve the (unphysical) x and y dependence of the thermal 
Lagrangian, then perform the differentiation and, subse- 
quently, set E = 0. In addition, the vanishing of the elec- 
tric field is required for the vacuum in order to remain 
passive according to assumption 4) . The appropriate ex- 
pression is simply obtained by replacing 

esB cot esB (37a) 

by the gauge and Lorentz invariant terms 

(es) 2 |y| coth^es(v/a; 2 +?/ 2 +x) 2 ^J cot (es (y/ x 2 + y 2 - x) 2 ^ 

(37b) 

in analogy to eq.(p8|). Executing the differentiation and 
setting E — 0, the thermal Q-factor reads 
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AQ(B,T) 



loo 

a 1 fds 



ttB 2 J s 





zsB coth esB — 1 esB 



sinh es-B 

oo 



■ cothesS 



(38) 



An analytic evaluation of eq. ( p8[ ) is only possible for cer- 
tain limiting cases. First, we consider a purely tempera- 
ture-modified vacuum with vanishing field strength. The 
proper-time integration then yields 

99 2 00 2 

A ^ B =°' T ) = S^E(- 1 ) n (^) 

n=l 

(39) 

where denotes a modified Bessel function and reflects 
the 5* 1 x M 3 topology of the finite-temperature coordinate 
space. 

The low-temperature limit is easily obtained, since the 
Bessel function is exponentially damped for large values 
of its argument. Thus, only the first term of the sum in 
eq.d39j) is important and we find 



AQ{B = 0,T 



_ 22 a 2 „ , ... 



(40) 



Here, we find the substantiation of why the Q-factor of 
the low-temperature velocity shift is simply given by its 
constant value at the origin in field space Q(B = 0, T = 
0) = C\ + C2 cq. ( p6| ) . Finite-temperature corrections van- 
ish exponentially in the low-temperature limit. 

The same conclusion holds for the parallel-plate 
Casimir configuration, the so-called Scharnhorst effect, 
due to the similarities of both effects concerning periodic 
boundary conditions: the corrections to the "universal" 
constant of eq.(^|) regarding this Casimir vacuum vanish 
exponentially with the plate separation a: T cx -. In this 
sense, the derivation of the Scharnhorst velocity shift is 
a trivial application of our light cone condition. 

Nevertheless, Q ceases to remain constant if we move 
noticeably away from the origin in field/parameter space. 
The high-temperature limit of the thermal correction 
eq.(|3|) to the Q- factor serves as an example. In this 
limit, the complete infinite sum has to be taken into ac- 
count. With some effort, one obtains 



AQ(T > m) 



22 

45^ 



k\ m 4 

Ir 1 



o 



to 



(41) 



where ki = 0.123 749 077 470... =const.. We immedi- 
ately find the complete Q-factor 



Q(T > m) = Q(T = 0) + AQ(T > to) 



11 a 2 



(m 2 



(42) 



which exhibits a rapid decrease cx 1/T 4 . Numerical anal- 
yses show that eq.(fi2|) is already valid below the e + e~ 
threshold where T/m < 1. 

With the aid of the finite-temperature VEV of the 
energy-momentum tensor 



90 



(N B + lN F ) T 4 diag(3,l,I,I) 



(43) 



where the integer variables Nb and N-p denote the num- 
ber of bosonic and fcrmionic degrees of freedom at a given 
temperature, we recover the well-known low-temperature 
result [i 



f 



44vr 2 
2025 ' 



T 4 



(44) 



where we used Nb = 2, iV F = 0, and eqs.@@||). In 
the high-temperature limit, the velocity of soft photons 
moving in a photon and ultra-relativistic e + e~ gas, i.e., 
./Vb = 2 and N-p = 4, is found by employing eq.([42|) 

121 



v = 1 - 



8100 
1 - 9.72.. • 10~ 7 



o 



— const. + O 



(&) ■ (45) 



Starting with an increase proportional to T for low tem- 
perature, the velocity shift approaches a constant value 
in the high-temperature limit. The magnitude of this 
constant indicates its subdominant importance even in 
extremely hot surroundings (early universe). 



D. Finite Temperature and Magnetic Fields 

A non-trivial interplay of the various effects is found in 
the domain of strong fields in hot surroundings. In any 
other limiting case, thermal phenomena decouple from 
the electromagnetic ones. Therefore, we evaluate the Q- 
factor eq.(Bq) in the limit B ^> B CI 



l 1 r . 



n=l 



(46) 



Subsequently, we take the limit T/m 3> 1 and find 

» ~ /™ \ la I B la, e m ( m 2 \ 
= -6n^. + l-^BT 2 + ° M 

(47) 

where k 2 = 0.213139199 408... =const.. Again, the 
first term in this expansion exactly cancels the zero- 
temperature strong-field contribution to Q found in 
cq.(|3l]) and we obtain for the complete Q-factor 
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2 a 2 I ( m 2 \ 

Q(T»m,B»B CI ) = -k 2 -, m+ 0^, (48) 

where we have introduced the convenient dimensionless 
variables B = -M- = anc j j 1 = H which satisfy 

B,f > 1. Note that the energy density (T 00 ) consists 
of three contributions: a thermal (eq.([43|)), a magnetic 
(\B 2 ) and a mixed part that can be obtained from the 
high-temperature limit of eq.(|36|). 

Finally, we arrive at the polarization and propagation 
direction averaged light velocity for a strong B-field and 
high-temperature modified vacuum 

v = l kiO?— -a 2 . (49) 

135 B 18tt t 2 27 K ' 

= 1 - 9.13.. ■ lCT 6 — - 2.75.. • 1CT 5 ^- - 4.21.. ■ 1CT 7 
B T 2 

One finds a minimal velocity shift at T 2 /B ~ 1.74 where 
\6v\ ~ 3.20 • 10 -5 . This magnitude corresponds to the 
velocity shifts we typically found in the present work for 
strong fields consistent with the one-loop approximation, 
e.g., in eq.(||). 

IV. CONCLUSIONS 

In this letter, we investigated light propagation in non- 
trivial QED vacua in the geometric optics approxima- 
tion. For a certain class of vacua fulfilling moderate as- 
sumptions, we derived the light cone condition averaged 
over polarization states which represents a generalization 
of the "unified formula" found by Latorre, Pascual and 
Tarrach ||. The latter was identified as the low-energy 
limit of our light cone condition. Especially, the "uni- 
versal constant" of eq.(|2|) turned out to be a derivative 
combination of the Lagrangian evaluated at the origin in 
field/parameter space. 

Within this conceptual framework, we calculated the 
velocity shifts induced by various vacuum modifications. 
In the low-energy limit, we recovered the well-known re- 
sults that were already perfectly described by the unified 
formula. As a high-energy example, we reproduced the 
findings of Tsai and Erber jl(| for intense magnetic fields 
using our comparably simple formalism. 

In the high-temperature limit, we observed that the ve- 
locity shift approaches a constant but comparably small 
value in contrast to the rapid increase cx T 4 for low tem- 
perature. Of course, this result is of minor importance 
from an experimental viewpoint. However, on the one 
hand, it at least confirms the consistency of the formal- 
ism, and on the other hand, as a non-trivial prediction 
of QED, it also indicates that the theory perfectly con- 
trols the velocity shifts. In the latter sense, this result 
confirms the consistency of QED. 

The only formally unbounded velocity shifts were dis- 
covered in the limit of strong fields at zero as well as high 



temperature. With regard to the validity of the loop ex- 
pansion, we expect higher-order or even non-perturbative 
Lagrangians to cure this problem. Then the light cone 
condition might serve as an indicator of consistency. 

Since the calculated velocities are phase and group ve- 
locities of soft photons, i.e., in the zero-frequency limit, 
the whole formalism is not immediately appropriate for 
discussing questions of causality that arise from the pos- 
sibility of positive velocity shifts, (v > 1 for negative 
energy densities, e.g., Casimir vacua); the maximal sig- 
nal velocity, namely, is evaluated in the infinite frequency 
limit. For indirect predictions, the reader is referred to 
refs. [p[JT5|] . The question of causality is extensively dis- 
cussed in ref. p[ . In concordance with these authors, 
let us just say that we find no grounds for a violation of 
the causal structure of space time by quantum vacuum 
effects. 
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